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There are hidden observables for inflation, such as features localized in position space, which do 
not manifest themselves when only one inflation trajectory is considered. To address this issue, we 
investigate inflation dynamics in a landscape mimicked by a random potential. We calculate the 
probability for bifurcation of the inflation trajectory in multi-stream inflation. Depending on the 
shape of the random bumps and the distance between bumps in the potential, there is a phase tran- 
sition: on one side of the critical curve in parameter space isocurvature fluctuation are exponentially 
amplified and bifurcation becomes very probable. On the other side bifurcation is dominated by a 
random walk where bifurcations are less likely to happen. 

PACS numbers: 



I. INTRODUCTION 

During multi-stream inflation (T]-|4] , the classical inflation trajectory bifurcates into multiple branches. Multi-stream 
inflation provides a possible explanation for features and asymmetries in the CMB, e.g. non-Gaussianities [U[3], cold 
spots in the CMB and great voids in the large scale structure [3J , and also have a web implication for eternal inflation 

m- 

Although multi-stream inflation is not necessarily embedded into string theory, the motivation is best understood 
in terms of the string landscape [SHE]. According to the string landscape paradigm, string theory is unique while 
the vacuum structure of string theory is extremely complicated. Based on semi-classical investigations of string 
compactifications one estimates that there are of order 10 100 or more meta-stable vacua. This forms the string 
landscape. The implications of the string landscape for inflation should be considered seriously. 

Over the past decade there has been a lot of work on inflation in the context of string compactifications (see e.g. 
[SHU] for some early papers and e.g [T2] for some recent review articles). A lot of the work has focused on specific 
string theoretical constructions which might yield inflation, and has not taken the landscape features of string theory 
into account. However, there have been a number of inflation scenarios motivated specifically by the string landscape. 
An incomplete list includes inflation in a random potential |13H15j , chain inflation |16H20| , particle production during 
inflation [21 -24], extra symmetry point encounter [251126) . meandering inflation |27j . multi-stream inflation and quasi- 
single field inflation H5] . 

In this paper, we focus on multi-stream inflation. It was argued that in the string landscape, where the inflation 
potential is extremely complicated, bifurcations in the inflation trajectory become possible or even frequent. But no 
quantitative justifications for this statement were previously provided. The aim of the present paper is to fill this gap, 
and bring a connection between the landscape and observation. 

The observational consequences fall into two possibilities: When the bifurcation probability is rather small, there 
may be only a few bubbles in the sky following exotic trajectories. In this case multi-stream inflation provides a 
possible explanation of the cold spot on the CMB, or a void in the large scale structure. On the other hand, when 
the bifurcation probability is large, extra fluctuations are not under control and thus such random potentials are 
not consistent with observations - except when they are extremely finely tuned 1 , or when the time of bifurcation is 
pushed to before the observable period of inflation (the final period of inflation during which scales which are currently 
observed exit the Hubble radius). 

The inconsistency for high probability bifurcations in a random potential comes in at least two aspects: Most 
importantly, different trajectories typically have different c-folding numbers, the e-folding number difference between 
different trajectories results in an extra fluctuation in the scalar type perturbations. Such a fluctuation with amplitude 
much greater than 10~ 5 is already ruled out on the scales of CMB. Also, if the different trajectories do not combine 
or reheat into the same radiation, there will be domain walls in between them, which causes a problem. 

In the string landscape the moduli space of low energy modes may have a very large dimension. For simplicity, we 
assume in this paper that only two fields participate in the dynamics during the inflationary phase. All other flat 
directions are taken to be stabilized. They do not evolve dynamically but could provide bumps in the effective two 



Here we are having in mind a random potential. If the bifurcation potential is not random (say, originate from spontaneous breaking of 
an approximate symmetry), common bifurcations do not necessarily require a fine tuning. 
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dimensional landscape. If there were more flat directions taking place in the inflationary dynamics, the bifurcation 
probability should be larger. Thus, in our present analysis we are actually under-estimating what the bifurcation 
probability would be in a more general case, perhaps by a huge factor. 

The paper is organized as follows: In Section 2, we model the landscape and calculate the bifurcation probability 
in the landscape analytically. Two cases are considered, either when bifurcations are dominated by an exponential 
growth of isocurvature fluctuation, or when they are given by a random walk of the isocurvature field. In Section 
3, we use stochastic equations to numerically verify the analytical approximations made in the previous section. We 
discuss our results in Section 4. 



II. BIFURCATION PROBABILITIES IN A LANDSCAPE 

In this section, we calculate bifurcation probabilities in the landscape. We show that there is a phase transition 
controlled by the height and the ellipticity of the bumps in the potential. 

To model the random potential, we focus on a double-field model. The fields are labeled ipi and <^ 2 - There is a 
overall slope in the (pi direction, such that if we remove the randomness of the potential, tpi becomes the inflation 
direction and ip2 has flat potential. To add random features to the potential, we introduce a perturbation U of 
the potential which is characterized by randomly located and separated bumps. These bumps are characterized by 
field values A p ipi and A p t/? 2 which denote the mean size and separation of the bumps in the ipi and Lp 2 directions 
respectively. An ellipticity parameter £ = A p ipi/ A p (f2 is defined to model the shape of the bumps. The amplitude of 
the bumps is taken to be the limiting value such that the slow roll condition in the inflaton direction is maintained. 

The equations of motion for the fields take the form 

(p x + 3H<p! + diVfa) + diV{<pi, <p a ) = , (1) 
^a+3Hip 2 + d i U(ip 1 ,ip2)=0, (2) 

where V(ipi) is the slow roll potential, subject to slow roll restrictions, and U(ipi, ipz) is a random potential. Shorthand 
notations d\ and di are used to denote the partial derivatives with respect to the two fields. We use a parameter A 
to denote the ratio 

\=^mjf)/\diV\. (3) 
When A -C 1, the slow roll condition in the inflation direction remains satisfied. 

A. Definition of a bifurcation 

One subtlety is the question of quantifying what is a bifurcation. We adopt the following two definitions: a 
bifurcation happens when two sample inflation trajectories (with the same initial conditions): 

1. End on different sides of a bump. 

2. Are much farther apart at the end of the inflationary phase than would be expected from quantum fluctuations 
alone. 

Numerically, we find that the bifurcation probabilities obtained from these two definitions are practically equal in the 
interesting region of parameter space. The two definitions give different results when the size of the bumps (both in 
terms of field extent and height) becomes too small to have any impact on the dynamics. In that case the probability 
drops to zero for definition [2] while can stay high for definition 1. As we will see, the size when the bumps in the 
potential become negligible is when they become smaller than the expected quantum fluctuations of the fields. With 
that insight we define bifurcations based on Definition [2] 2 . 



2 Definitions 1 and 2 also differ when the field extent of the bump is as large as the field distance which the inflaton field moves in of 
order of a few e-foldings. However, in this case we cannot speak of a random potential any more inasfar as the dynamics relevant to 
cosmological observations is to be considered. Thus, here we are not interested in this situation. 
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Bifurcation from instabilities 



Consider two trajectories of ip%, denoted by ip 2 



(•-i) 



(A) 
¥>2 



and 



Let the deviation of these two trajectories as 5 = 



ip 2 ' ■ Then as long as 8 is smaller than A p Lp 2 , one can expand U such that S satisfies 

5 + 3H8 + (d%U)5 = . 



(4) 



Now we discuss general features of the potential U. It follows from the classical equation of motion and from the 
normalization of the density fluctuations produced during inflation that during one e-folding of inflation (fx direction 
rolls a distance of order 10 5 H. 3 This distance is extremely long from the point of view of the field (p 2 , because of 
the smallness and randomness of the potential U. Thus it is a good approximation to describe the dynamics of ifi2 as 
taking place in a time-dependent potential given by 



U((pi,(p 2 ) = U((fi(0) + (pit,(p 2 ) , 



(5) 



where <fi(0) is the initial value of the field and (fx can be treated as a constant. 

Over a time interval of a few Hubble times, the field trajectory can pass through a lot of bumps. The number of 
bumps the field trajectory crosses during a time interval t is about tfxt/ A p <fx- Thus, on time scales of order A p (fx/<P, 
the d 2 parameter in equation Q will typically change its sign. To model this behavior, we can use the following 
approximation for the form of the second derivative of the potential U((f2,t) from jHJ: 



■ sin 



x -i x I (6) 

A p tp 2 \AptpiJ 

When the sine function change its sign, the equation Q changes its behavior from oscillation to exponentially growth. 
Now the question is whether the overall behavior shows exponential growth or not. Qualitatively, we see that when 
the sine function changes its sign slowly enough, there is enough time for growth to take place. On the other hand, 
when the sign function change its sign extremely quickly, the equation Q can be treated adiabatically thus there 
should be no growth. 

To be a bit more precise, if we neglect the Hubble friction term for a moment (we will return to the full treatment 
at the end of this subsection) , the equation Q becomes 



S + asia(fit)5 = Q 



(7) 



and thus has the form of the Mathieu equation, the same equation which describes preheating after inflation (see 
[301 [31] and [32] for a recent review). It will have exponentially growing solutions when 



\a\//3 2 > 0.45 



(8) 



(as shown in figure [I]) , and non-growing solution otherwise. Note that ^ is the condition for the time scale of the 
change of the sign of the frequency to be longer than the intrinsic time scale of oscillation. 

Inserting the parameters describing the size of the random potential, the condition ^ to have growing solution in 
the <p 2 direction becomes 



A p ip 2 



> 0.45 



2nipx 

Aplf! 



(9) 



i li 

Note that the slow roll condition in the inflaton direction implies |c?iV| — |3if^i|. Making use of the result = 
H 2 / (2ir(p) for the power spectrum of cosmological perturbations £ to eliminate ip\, the above inequality becomes 



A p ipx 
H 



> 



1/2 
C 



(10) 



When this inequality is satisfied, the solution turns out to be exponentially growing. In order for the exponential 
growth rate to be large on compared to the Hubble expansion rate we need 



\£dxV 



A p ip 2 



> H 2 



(11) 



This is because, the rolling distance of tp\ per Hubble time is ip/H ~ (H 2 /ip) 1 X H ~ P, 
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FIG. 1: The behavior of equation 5 + asm(/3t)8 = 0. Here /3 = 1 is chosen without loss of generality because varying ft is a 
rescaling of time. The left panel has a//3 2 = 0.45, and the right panel has a/0 2 = 0.46. A sharp transition is observed: for 
a = 0.45, the amplitude of 8 does not grow, while for a — 0.46, the amplitude grows exponentially. 



This requirement can be rewritten as 



Thus when 



H 2ttP c 1/2 



<=*<^. (13) 



a perturbation in the isocurvature direction will grow until its size reaches the value A p tp 2 when bifurcations become 
common. Note that in terms of £, both the inequalities give lower bounds. 

Now consider the full equation (4| including the Hubble friction term. Experience from studies of preheating (see 
[32] for a recent review) teach us that an instability derived in the absence of Hubble friction persists when restoring 
the friction term: via a field rescaling with a power of the cosmological scale factor the friction term can be eliminated, 
and the equation of motion of the rescaled variable remains of the type described by Floquet theory and showing an 
instability. 

If there is to be an observable bifurcation, we require that after about 10 e-foldings the growth factor times H/ (2ir) 
(which is the typical magnitude of quantum fluctuations in a de Sitter phase and which we take to be the initial value 
of 5) becomes greater than the separation A q ip 2 between bumps in the random potential U . This condition can be 
analyzed numerically, and the result is shown in Figure [2] and Figure [3] as a lower bound on £ to get bifurcation 
leading to amplification of the fluctuations. In both figures the area above the curve yields the parameter space for 
which bifurcations occur and have an observable effect. 

From Figure [2j we observe that for small A p ipi/H, one can safely neglect 3Hip 2 . This is because in that case the 
oscillations have a much higher frequency than Hubble friction. On the other hand, when A p ipi/H is large, neglecting 
3H(p2 results in an over-estimation of £ by a factor of 1.7. Figure [3] compares the lower bound for bifurcations obtained 
using different values of A. 



C. Correction factors 



In the previous subsection, we discussed the basic picture of amplified bifurcations in a random potential. However, 
the model we used is over-simplified, in the sense that there are factors of order one (or order 10) which were ignored. 
Since we want to relate our analytical analysis to numerical simulations we must discuss the various correction factors 

First, as will be discussed in Section [ill A[ the numerical simulations are performed with the maximal amplitude A 
of the perturbation U of the potential being a function of the bump size, i.e. A(A p <^i , A p ip%) (chosen such that the 
motion in the ipi direction is at the borderline of slow rolling). In the analytic calculation this corresponds to working 
with A = 1/(100- 

Second, there is a correction factor along a single inflation trajectory. We have used A p ipi and A p (p 2 to denote 
the average distance between bumps in the (pi and ip 2 directions. However, the criterion for amplified bifurcation is 
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FIG. 2: Comparison between the results neglecting the 3Hip2 term (lower curve) and the full result (upper curve). Here the 
value A = 0.1 is used. The estimate neglecting friction overestimates the bifurcation probability by a factor of up to 1.7 in the 
variable £. 
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FIG. 3: The lower bound on £ in order to have bifurcation of the amplification type. Values of £ above the red, green, blue lines 
(top line is the red line, bottom line is the blue line) correspond to regions in parameter space where exponential amplification 
occurs for values of A =0.02, 0.10, 0.50 respectively. In these regions of amplification, the bifurcation probability becomes of 
order one. 



actually not controlled by the average distance, but rather the longest distance between bumps of the same sign in 
the potential. 

To put it more explicitly, along an inflation trajectory, we have assumed that the derivative of the potential in the 
isocurvature direction d^U changes its sign when the inflaton rolls an average distance A. p if±. However, a random 
potential is not exactly periodic. As a result, on a distance of order nA p tpi (where n is a positive integer), the sign 
of &2U as a function of time will be a sequence with length n, like (+, +, — , +, — , — , —,—,+,••• ). The probability 
of bifurcation is determined by the longest sub-sequence without a sign change (4 in the above example). The 
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determination of this critical distance is a well defined mathematics problem. Here we do not give an analytical 
solution. However, numerically, we find that the expectation value (as an ensemble average) of the longest sub- 
sequence without sign change (n su b) can be approximated by 



{n su b) 



1.51og(n) 



(n > 1) 



(14) 



If we are interested in AiV e-foldings of inflation, n can be written as 

ipAN AN H 



(15) 



For example, if we take AN = 10, use the CMB normalization of the curvature power spectrum, and take A p ipi/H = 
100, we have n = 324 and thus (n su b) = 8.7. For A p tpi/H = 10 3 , we have (n su b) = 5.2. 

This means that if A p tpi is the mean size of the bumps in a random potential, then we should compare our numerical 
simulations to analytic results with bump taken to have size (n su b)A p ipi. Alternatively, if A p ipi is the size of the 
bumps in the analytic analysis, we must then compare the analytical results to simulations in a random potential 
with bumps of mean size l/(n su b) A p (p\. This (n S ut>) factor yields a correction of the parameters A p tpi, £ of the form: 



Apipx 
_ A p ipi 



i= 



<n au b> A p <,02 ' 



(16) 
(17) 



and should be used in equations ([6| ~ (13). The bifurcation probability using these first two correction factors is 
plotted in Figure |4j For large (n su b), we have used the above approximated formula (14), and for small (n su b), we 
have used the explicit result. 
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FIG. 4: Probability of bifurcation from analytical calculation simulations. We have rescaled data according to equation 
Also, a color plot is made for the purpose of comparison with numerical results, with the same set of parameters. 
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Third, in the previous subsection we were comparing only two trajectories. However, there are exponentially many 
trajectories (local Hubble volumes) created during inflation. Thus there may be rare bifurcations in some exceptional 
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trajectories. To take those exceptional trajectories into account, note that the quantum fluctuations are nearly 
Gaussian and thus given by the following probability distribution 

Prob oc exp(~Atp 2 /2a 2 ) , (18) 

where a is the variance of the fluctuations. On the other hand, there are exp(3A7V) trajectories, where AN is the 
c-folding number difference between the bifurcation point and the start of observable inflation. Thus, typically, the 
largest fluctuation among these exp(3A7V) trajectories is A(p ~ aVSAN. For example, if we consider AN = 10, 
where bifurcation takes place on the smallest observable scales on CMB, the rare bifurcations could take place even 
when Aptpi and A p tf2 are 5 times greater than their critical values. 

As mentioned in the introduction, there are big observational differences between the case of rare bifurcations and 
that of common bifurcations. The case of rare bifurcations could be consistent with current CMB observations and 
could in fact be responsible for a position space non-Gaussianity. However, if bifurcations become common, fine 
tuning is needed to make the scenario consistent with observations. This is because in a random potential, different 
trajectories commonly lead to different number of e-folds of inflation. The difference in e-folding number is translated 
into temperature anisotropics in the CMB via the SN formalism, which will not be consistent with observations when 
SN > 10" 5 . 



D. Bifurcation from a Random Walk 



Now let us proceed to consider the case in which bifurcations do not lead to an amplification of the isocurvature 
fluctuations. In this case, we define At as the time duration for the field to cross the bump in ipi direction. Thus 

At = Apipi/ip. (19) 

During the time interval At, the quantum fluctuation in the if2 direction is 

A q ip 2 = ^^HAt . (20) 

Z7T 

The bifurcation probability during this At time is determined by the ratio 

Prob(iJAt) = (21) 

Ap(^ 2 

when A q tf2 < A p (f2- Otherwise, when A q ip2 > A p <p 2 , bifurcations happen frequently, but on the other hand the 
bumps are too small to lead to interesting consequences. 

Now, consider a world line of a comoving observer. The probability for the observer to see bifurcation(s) during N e 
e-foldings is 

N B /(HAt) 

Prob(7V e ) = t 1 ~ VnProb(HAt)} . (22) 

n=l 

The y/n factor comes from the random walk accumulation. To have a considerable bifurcation probability, we need 

3/2 

Prob(ff At) > 1 . (23) 
Inserting numbers, we obtain the following condition for bifurcation 

Apipx 



2 / N e V 

3 \HAt J 



H 



< 



at 3 / 2 /: 

iVe 4 (24) 



\ 6tt 2 P c 



1/2 



For 10 e-foldings, which corresponds to observable scales on the CMB, and for £ — 1, the above inequality is saturated 
for the value A p ipi/H ~ 100. 

In this subsection, we have not considered the classical dynamics of the ip2 field. Considering the ip 2 field may 
encounter a bump and bounce back, we may have over-estimated the bifurcation probability. Moreover, if we impose 
the slow-roll constraint e -C 1, bumps formed of Lagrange polynomials will have a negligible effect on the fields 
compared to quantum fluctuations when A q ip 2 / A p ip2 ~ 1. Hence we do not expect any bifurcations if we stick with 
the definition [2] The conclusion of the current subsection is that, without amplification, bifurcation is less probable 
to happen in a way which has observational consequences. 



III. SIMULATION FOR STOCHASTIC INFLATION 



Numerical simulations were performed to obtain the probability of bifurcation in random potentials. The dynamics 
of our fields are created using Starobinsky's stochastic approach [33] to take into account the quantum fluctuations. 
In this approach, the equations of motions for tpi and (p 2 take the form, 

<h(t) + 3H<pi(t) + V m = ^H b ' 2 m {t), (25) 



(p 2 {t) + 3^ 2 (<) + V V2 = ^H 5 / 2 m (t). (26) 

Here the rji terms are stochastic sources which obey independent Gaussian distributions and are normalized to give 

(r H (t)r h (t'))=6(t-t'). (27) 
To evolve the fields numerically we discretize time in intervals At = t n — t n -\ so that the delta function becomes 

<f(*-*0">^J- (28) 

We do not use the approximations \(p\\ <C |3_ff<£>i| during the simulation but we did make use of the approximation 
that the energy density is dominated by the potential in order to write H — 



A. Random Potential 



Our random potentials can be written as, 

Vteii + Ubuw) , (29) 

where V is a background potential that has a small tilt in the <p\ direction. We chose V — \m 1 Lp\ for the numerical 
work. The second term U is a set of random perturbations which are created by discretizing field space (ipi, f 2 ) into 
a lattice with separation 4 A p ipi and A p (p 2 and assigning a random number from the interval [—A, A] at each point, 
and then performing a two dimensional interpolation to obtain a continuous function. 

The value of A(A p (pi, A p ip 2 ) > represents the amplitude of the perturbations and depends on the lattice spacing 
because we require that the full potential V + U does not lead to violation of the slow-roll condition on the e parameter. 
This constraint can be made explicit when the form of the perturbations are known. In our case the perturbations are 
interpolations between the random values associated to each lattice point using a second order Lagrange polynomials. 



B. Results 



The probability of bifurcation is expressed as a two dimensional plot. The ellipticity £ and Aipi/H are used to 
quantify our parameter space. 

To achieve a better understanding, two types of potentials are considered: a toy potential with regularly arranged 
bumps (the left panel of figure [5]) , and a completely random potential (the right panel of figure [5]) . The toy potential 
is faced with the same problem as the analytic model if we want to compare it to the random potential: we must 
take into account the factor (n su b) and rescale the parameter axes accordingly. The latter is what we want eventually. 
Note that the figures are only for illustrative purpose, and the actual size of the bumps are much smaller. The results 
are shown in Figure [6] The plots can be compared with Figure |4j which is obtained from the analytic model. 

The shape of the maps from our numerical results and based on the analytical analysis agree with each other 
when Aptpx/H < 10 4 . When A p (p±/H > 10 4 , then in the numerical simulations the bifurcation probability decreases, 
whereas the analytical calculation does not have lead to such a decrease. The difference arises because in the analytical 



4 A given lattice spacing will represent a potential which is dominated by perturbations whose mean size is of the lattice spacing, hence 
the use of A p <fii as notation. 
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FIG. 5: An illustration of potentials with regularly placed bumps (left panel) and random bumps (right panel). The size of 
bumps in the actual calculations are much smaller. 

analysis we are assuming that the trajectory encounters a large number of bumps. However, when A p <pi/H > 10 4 , the 
above assumption becomes invalid. We also note that the analytical result under-estimates the bifurcation probability 
by a factor of about 10. 

Using the numerical code, we can also study explicitly what observational problems occur in the case that the 
bifurcation probability is large. In Figure[7j we show two runs without (left panel) and with (right panel) bifurcations. 
It is shown that when the bifurcation rate is large, the distribution of e- folding numbers in local patches of the universe 
becomes highly non-Gaussian, which is not consistent with current observations. 

IV. CONCLUSION AND DISCUSSION 

To conclude, we have considered the probability of bifurcation of the inflaton trajectory in a (possibly stringy) 
inflationary landscape. A phase transition from bifurcating dynamics to absence of bifurcations is identified. The 
phase transition line has been found as a function of the mean distance of the bumps (A p (pi/H) and of the ellipticity 
parameter A p <pi/ A p ip2- 

There are a number of interesting questions in this direction which we have not addressed in the present work. For 
example, the analytical analysis has under-estimated the bifurcation probability by a factor of about 10. Perhaps we 
are simply using too strong a criterium for bifurcation. But it is also possible that there are other effects leading 
to bifurcation, that we do not understand currently. Also, in the analysis leading to the results of Figure [4] we only 
considered the effect from exponential growth of isocurvature perturbations. The random walk may also contribute 
in the region A p ipi/H < 10 3 . A combined analysis is needed to get a more precise analytical prediction. 

As another example, as we mentioned in the introduction, a string landscape could have a large number of modular 
fields - maybe 0(100). A lot of them may be light and have dynamics. We only considered two of them. It remains 
interesting to study the behavior in a many-dimensional random potential, and to determine how the bifurcation 
probability scales as a function of dimension of field space. 

Also, a many-dimensional field space opens up the possibility of bifurcation in isocurvature directions. The bifurca- 
tion in isocurvature directions could have completely different predictions compared with bifurcations in the curvature 
direction, and this issue has not been investigated. 

In the present work, we have used the separate universe approximation and neglected the field gradient terms. It 
would be interesting to add those terms back and do a full simulation in position space. Our current analysis does not 
consider the effect from the boundary regions between bubbles following different trajectories. These regions should 
become domain walls for a period of time. In the case of rare bifurcations, the trajectories could recombine, and thus 
there is no domain wall problem. However, it remains interesting to see whether these wall-like objects could lead to 
observationally interesting predictions, or rule out the whole model. 

Finally, in the case of rare bifurcations, i.e. in the parameter regime on the boundary of the black regime in Figures 
|4]and[6j it would be interesting to see what kinds of bifurcations are more natural in a landscape. The investigation 
of this question may lead to a more precise correspondence between the string landscape and observations such as the 
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Toy potential 




FIG. 6: Probability of bifurcation from numerical simulations. The upper figure is the original result of the toy potential with 
regularly spaced bumps and the bottom left panel is the same result taking into account the correction factor from (n su b). The 
bottom right panel is the result for a random potential. 
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FIG. 7: Probability distribution for the fluctuations based on the 5N formalism. The left panel is for trajectories without 
bifurcation, and the right panel is for trajectories with a large probability of bifurcation, without symmetry protection or fine 
tuning. 
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CMB cold spot. 
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